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Abstract —In this paper, we study the simultaneous wireless 
information and power transfer (SWIFT) in downlink multinser 
orthogonal freqnency-divlslon multiple access (OFDMA) systems, 
where each user applies power splitting scheme to coordinate the 
energy harvesting and secrecy information decoding processes. 
Assuming equal power allocation across suhcarriers, we formu¬ 
late the problem to maximize the aggregate harvested power of all 
users while satisfying secrecy rate requirement of individual user 
hy joint suhcarrier allocation and optimal power splitting ratio 
selection. Due to the NP-hardness of the problem, we propose 
two suboptlmal algorithms to solve the problem in the dual 
domain. The first one is an iterative algorithm that optimizes 
subcarrier allocation and power splitting in an alternating way. 
The second algorithm is based on a two-step approach that solves 
the subcarrier allocation and power splitting sequentially. The 
numerical results show that the proposed methods outperform 
conventional methods. It is also shown that the iterative algorithm 
performs close to the upper bound and the step-wise algorithm 
provides good tradeoffs between performance and complexity. 

Index Terms —Physical-layer security, simultaneous wireless 
information and power transfer (SWIPT), energy harvesting, 
orthogonal frequency-division multiple access (OFDMA). 

I. Introduction 

Orthogonal frequency division multiplexing access 
(OFDMA) gains its popularity and has become a leading 
multiple access candidate scheme for beyond fourth/fifth 
generation (4G/5G) wireless systems, due to its flexibility 
in resource allocation and robustness against multipath 
fading. It enables efficient transmission of various data traffic 
by optimizing power, subcarrier, and bit allocation among 
different users. 

Due to the broadcasting nature of wireless channels, secu¬ 
rity is a crucial issue in designing wireless communication 
systems. As a traditional method, cryptography encryption 
dominates the upper layers mainly by increasing the complex¬ 
ity has been introduced in every layer but physical layer in 
the standard five-layered protocol stack. Thus physical-layer 
security is an important complement to the other security 
approaches. 

A great deal of studies have been devoted to the 
information-theoretic physical-layer security |[T|-||9l. For ex¬ 
ample, in a, 0 , resource allocation for physical-layer se¬ 
curity considerations was studied for multicarrier systems. In 
0, artificial noise was considered for physical-layer security. 
However, the artificial noise based methods for physical-layer 
security mainly lie on the spatial degrees of freedom offered 
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by multiple antennas to degrade the channel of the eaves¬ 
droppers. Alternatively, simultaneous wireless information and 
power transfer (SWIPT) also becomes an important solution 
to improve the energy utilization for wireless networks by 
prolonging the lifetime of wireless nodes and draws a great 
deal of research interests m-ini. Energy harvesting wireless 
networks are potentially able to gain energy from the wireless 
environments. The prior work Qol studied the performance 
of SWIPT in the receiver that can decode information and 
harvest energy for the same received signal, which may be 
not realizable however. Two practical schemes, so-called time 
switching and power splitting, were proposed in ifTTIl . ifT^ as 
practical designs. With time switching applied at a receiver, 
the received signal is either processed for energy harvesting or 
for information decoding. When the power splitting is applied 
at the receiver, the received signal can be split into two streams 
with one stream processed by the energy receiver and the other 
processed by the information receiver. The authors studied 
the flat-fading channel variations in SWIPT in ca, where 
dynamic power splitting was applied in the systems. Two 
SWIPT schemes in OFDMA with different configurations and 
corresponding resource allocation problems were studied in 

na. 

A handful of works have studied SWIPT for physical- 
layer security, usually considering some receivers decode 
confidential information and the rest receivers harvest energy 
(also known as the separated receiver model) Ha-El. 
These works are mainly motivated by the dual use of the 
artificial noise, i.e., artificial noise is used to interfere with the 
eavesdropper for secrecy information receivers and acts as the 
source of energy harvesting for energy receivers. Such method 
is efficient for the separated receivers, however not for the 
co-located receivers where the receivers can simultaneously 
receive secrecy information and harvesting energy. Greatly 
different from the existing solutions for physical-layer security 
(such as artificial noise and beamforming), we consider a co¬ 
located SWIPT system by using power splitting scheme which 
is an “SWIPT” way against eavesdropping in an OFDMA 
system. Specifically, if subcarriers are preferable to transmit 
secrecy information, the user may split more received power 
for information decoding and, on the contrary, if subcarriers 
are easily eavesdropped, the user splits more received power 
for energy harvesting, which helps the systems to fully utilize 
both spectrum and energy of the easily eavesdropped subcar¬ 
riers that are traditionally difficult to utilize. 

One challenge is that power splitting should be per¬ 
formed before OFDM demodulation in practical OFDMA- 
based SWIPT systems. Thus each user should split the power 
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of all received information on all subcarriers at a same ratio 
instead of dynamic ratio. The complexity of designing power 
splitting ratio on all subcarriers with the same ratio is much 
greater than that of splitting on each subcarrier with dynamic 
ratio, because in the former case, the power splitting ratio of 
each user couples all subcarriers in the rate expression. 

In this study, we consider the secrecy-rate required downlink 
multiuser OFDMA networks, where all users apply power 
splitting scheme (of a same power splitting ratio at all subcarri¬ 
ers) to coordinate energy harvesting and information decoding 
processes. By assuming equal power allocation at subcarriers, 
our goal is to maximize the aggregate harvested power of all 
users while satisfying the secrecy rate constraint of each user 
by jointly optimizing the subcarrier allocation and designing 
the power splitting ratio. 

We formulate the problem with power splitting applied 
at each receiver for practical application (P-PA) as a mixed 
integer programming problem and NP-hard. Since the optimal 
solution is difficult to obtain, we introduce two suboptimal 
algorithms with polynomial time complexity. We first propose 
an efficient iterative algorithm to find the power splitting ratio 
and subcarrier allocation in an alternating way. To further 
reduce the complexity, we also propose a two-step algorithm 
that first obtains the optimal subcarrier allocation policy and 
then solves the optimal power splitting ratio. It is shown to 
tradeoff the complexity and performance. Numerical results 
show that the proposed iterative algorithms perform close to 
the performance upper bound and both proposed algorithms 
outperform the heuristic methods. 

The rest of this paper is organized as follows. In Section II, 
we formulate the problem. In Section III, we propose two so¬ 
lutions. We study the case of statistical CSI of eavesdropper in 
Section IV. In Section V, the performance of the two schemes 
are evaluated via numerical results. Finally, we conclude with 
a brief summary of our results in Section VI. 

II. System Model And Problems Formulation 

In this paper, we consider a downlink OFDMA network, 
which consists of one base station (BS) with one antenna, 
K mobile single-antenna users, over N subcarriers and one 
single-antenna eavesdropper attempting to wiretap informa¬ 
tion from all subcarriers. It is trivial to extend to the non- 
cooperative multi-eavesdropper scenario, since the overall 
eavesdropped rate is the maximum rate of the multiple non- 
cooperative eavesdroppers. Thus the proposed algorithms are 
also applicable if we select the best eavesdropper link among 
multiple eavesdroppers on each subcarrier, i.e., the eavesdrop¬ 
per with the highest decodable information rate. Each user 
communicates with BS and demands a secrecy rate that is no 
lower than a constant Ck > 0, for all 1 < /c < AT. Here we 
assume that equal power allocation is performed by the BS 
over all subcarriers for simplicity. This is reasonable since the 
gain brought by power adaption is limited in OFDMA systems 
(HSl-El. Each receiver is considered to split the received 
signal into two signal streams, with one stream to the energy 
receiver and the other one to information receiver. 

The considered OEDMA-based SWIFT method for 
physical-layer security can be applied in various scenarios. 


such as a home internet-of-things. In the considered example, 
the wireless devices such as phones and tablets are simultane¬ 
ously receiving confidential information and harvesting energy 
from the wireless access point (like Wi-Ei or femtocells). The 
wireless access point uses OEDMA to transmit signals (it 
is supportable in many standards). However, the neighbours 
in/around the building attempt to eavesdrop the secrecy infor¬ 
mation. 

We assume that all users are legitimate users and they have 
their own data transmission with the BS so that the BS can 
obtain full channel state information (CSI) of users. Let hk,n 
denote the channel gain of user k on subcarrier n, and /?„ 
denote the channel gain of the eavesdropper on subcarrier n. 
We also assume that each /3„ is independent and identically 
distributed (i.i.d.) Rayleigh fading channel. 

Let pn represent the fixed and equal power allocated on 
subcarrier n. The received signal at user k is processed by 
a power splitter, where we assume a ratio pk of power is 
split to energy receiver and the remaining 1 — of power is 
split into the information receiver for OEDM demodulation. 
We have 0 < < 1, Vfc. Note that power splitting is 

performed in analog domain before the digital domain where 
OEDM demodulation is processed. Thus, due to this hardware 
limitation, each user has to harvest the received signal with a 
same power splitting ratio on all subcarriers. 

With the full CSI of eavesdropper known to the BS, the 
achievable secrecy rate at subcarrier n of user k is given by 

HD 

= log, (^1 + ^ ^ 

( 1 ) 

where [•]+ = max{-,0}, cr^, rk,n and re,n are the power of 
additive white Gaussian noise, the achievable information rate 
of user k and the eavesdropper, respectively. Note that the 
full CSI of eavesdropper case is practically valid in following 
scenarios; (i) the eavesdropper is active in the network so 
that the BS can monitor its behavior and obtain its CSI; (ii) 
interestingly, as stated in || 23 , even an passive eavesdropper’s 
CSI can be obtained through its local oscillator power inadver¬ 
tently leaked from the receiver RE front end using the methods 
in lED, GS]; (iii) the legitimate users and the eavesdropper 
belong to different networks in today’s heterogeneous network, 
then the BS can coordinate with the eavesdropper’s serving 
network to obtain the CSI, since the eavesdropper is the 
legitimate user of different network or service. This is referred 
to as coordinated multi-point (CoMP) transmission in 3GPP 
LTE-A. The assumption is widely adopted in the physical- 
layer security literature (e.g., ED, Gbl-iEll). 

The secrecy rate of user k is given by 

N 

rk = '^Xk,nrln^ ( 2 ) 

n—1 

where we let Xk,n denote the binary subcarrier allocation 
variable, with Xk,n = 1 indicating that subcarrier n is assigned 
to user k and Xk,n = 0 otherwise. Note that if ^k,n = 0 






3 


for any subcarrier n, i.e., such subcarrier is not assigned to any 
user, then it is used to transmit power to users only. 

With the conversion efficiency of the energy harvesting 
process at each receiver denotes as 0 < c < 1, the harvested 
power of user k is thus given by 

N 

— CPk ^ ^ Pn^k.n ■ (3) 

n—1 


The goal of the considered problem is to find the optimal 
subcarrier allocation and power splitting ratio to maximize 
the total harvested power (for the purpose of uplink trans¬ 
mission for example) while satisfying the individual secrecy 
rate requirement for each user. This practical application 
optimization problem can thus be expressed as 


K N 


max 

{X.p} 

C ^ ^ Pk ^ ^ Pnhk,n 

k—1 n—1 

(4) 


K 


s.t. 

'^Xk^n < l,Vn 

(5) 


k^l 



^k,n ^ ,V/c,7Z 

(6) 


0 < Pk E 1, V/c 

(7) 


N 



J2^k,nrln >Ck,^k 

(8) 


n—1 


where X = {xk,n} and p = {pk}- The constraints in (|5]l and 
(l6]l enforce that each subcarrier can only be used by one user 
to avoid the multi-user interference. 


III. Proposed Algorithms 

The formulated (P-PA) is nonconvex due to the binary 
subcarrier variable Xk,n, finding the optimal solution is usually 
prohibitively due to the complexity. However, according to 
GUI, the duality gap becomes zero in multicarrier systems as 
the number of subcarriers goes to large and the time-sharing 
condition is satisfied. Thus the optimal solution of a nonconvex 
resource allocation problem in multicarrier systems can be 
obtained in the dual domain. 

Nevertheless, as we will discuss later, the traditional La- 
grangian decomposition cannot be directly employed to de¬ 
compose the problem into parallel subproblems with each 
subproblem corresponding to one subcarrier. This is because 
the power splitting ratio pk appears in the rate expression 
and couples the subcarrier assignment variables. As a result, 
solving (P-PA) is nontrivial though the dual method is used in 
this paper. In this section, we propose two efficient suboptimal 
algorithms. 


A. Iterative Algorithm 

We define T as all sets of possible X that satisfy Q and 
(|6]l, IZ as all sets of possible p that satisfy 0 < < 1. 


The Lagrangian function for (P-PA) is given by 
L{p,X,p) 


K 


N 


K / N 

= Pnhk,n + ^Pk ( X! “ ^k 

\n—l 


k—1 n—1 

K N 


/c=l 


K 


= EE i^CPkPn^k^n ^k.nPk'^k.n') -J2PkCk, (9) 


k—1n—1 




where p = [pi, p 2 ) Pk] are the Lagrange multipliers. The 
dual function is then defined as 


gip) = ^ max L{p, X, p). 

X^T 


( 10 ) 


The dual problem is thus given by min^ g{p)- For the max¬ 
imization problem in (fTOl i. the Lagrangian function cannot be 
decomposed into N subproblems, because the power splitting 
ratio Pk has to be computed considering all subcarriers that 
are assigned to user k, instead of one specific subcarrier. 

Thus, for given dual variables p, we can obtain a subop¬ 
timal solution by iteratively optimizing X with fixed p, and 
optimizing p with fixed X. The process is repeated until both 
X and p converge, which is known as the block coordinate 
descent (BCD) method lIUTI . 

To solve X with fixed p, suppose that subcarrier n is 
assigned to user k, we have 


N 


K 


L — ^ ^ Lfi ^ ^ PkCk'i 


k=l 


where 


K 


Lji — CPn ^ ^ Pk^kj, 


Pkrk,r. 


k=l 


Thus, the subproblem is given by 


max L„(p, JA„,/r) 

Xn^l 


( 11 ) 


( 12 ) 


(13) 


which can be solved independently. By maximizing each 
the optimal X can be obtained as 


^k,n 


1, k = k* = argmaxfe L„ 
0, otherwise. 


(14) 


To solve p with given X, the problem can be decomposed 
into K subproblems with each corresponding to one user 
since each pk is fixed in this process, which can be solved 
independently. The subproblem at user k is given by 

N 

max Lkipk) = {CPkPnhk,n + Xk,nPkrl „) , (15) 

i 

n—1 

and we have 

K K 

L = Y,Lk-Y,PkCk. (16) 

k=l k=l 

Applying the Karush-Kuhn-Tucker (KKT) conditions, we have 
each has to satisfy 


dLk 

dpk 


N 


= E 


CPn^k^n 


pkXk^nhk.nPn 


In 2 {hk,nPn (1 - Pfc) -I- 0-2) 


= 0 . 


( 17 ) 
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initialize p and p. 
repeat 
repeat 

Solve assignment variables X according to d and 
compute L according to (|9]l. 
for each user k do 

initialize ^ = 1 and = 0. 
repeat 


r,UB 




Algorithm 1 Proposed Iterative Algorithm for (P-PA) 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 
19 


Set Pk = h {Pk 
Compute according to (fT7] i. 

if > 0 then 

Set Pk = Pk- 

else 

Set p^^ = Pk- 

end if 


until 


dL^ 


< e, where £ is a very small constant 


dpk 

for controlling accuracy. 

end for 

until Lagrangian function converges. 

Update p by ([19]) according to the ellipsoid method, 
until p converge. 


Unfortunately, there is no closed-form expression for the 
optimal Pk. However, since Lk is a concave function of pk, and 
monotonically decreases as pk increases, we can adopt 
the bisection search method to solve p^ over 0 < pk < 1- 
Nevertheless, an asymptotic solution can be obtained by 
considering a high received signal-to-noise (SNR) scenario, 
i.e., tr^ —0. We have 


1 - Pfc = 


Pk 

Cln2 

^ l^n=l 


1 ^k,n 
Pnhk,n 


(18) 


In dH, is a constant in each iteration, '^n=iP^^k,n is 
user fc’s total received power which is also a constant, and 
'^n=i ^k,n is the number of subcarriers allocated to user k. 
Thus, we can conclude that 1 — pk, the ratio of the power 
splitting into user fc’s information receiver, is proportional to 
the number of subcarriers allocated to this user in high SNR 
scenario. 

With the fixed p, the optimal X* can be obtained by (O. 
The optimal value of the objective function can be increased by 
optimizing X via ( fTTl i. Then, with the fixed X*, the optimal 
p* can be obtained. Thus, the above process can be iterated 
until the optimal value of the objective function ceases to 
increase. 

Finally, according to 1^ . the dual function in ([TOl i is always 
convex. By simultaneously updating p, we can solve this 
problem by the subgradient method. The dual variables p are 
updated in parallel as 




(t) 

Pk 


N 


Ctk 


c„-Y, 


Xk,nrk.n 


n—1 


,Vfc. (19) 


The above iterative algorithm to solve (P-PA) is summarized 
in Algorithm 1. For this algorithm, the complexity mainly lies 


in step 19). As each pk is obtained individually by the bisec¬ 
tion search, the complexity of steps 6)-15) is 0{K). Hence, 
The complexity of steps 4)-16) is given by 0{K+KN). Next, 
the complexity of subgradient updates is polynomial in K 1^ . 
The overall complexity is given by where 

g is a constant and equal to 2 for the ellipsoid method. 

B. Step-Wise Algorithm 

Since the complexity of the above algorithm becomes 
unfavorable for practical application with the increase of K 
and N, we also propose a simpler suboptimal algorithm in 
this subsection. 

To begin with, we first formulate a problem by assuming 
that the power splitting can be designed differently on each 
subcarrier at each receiver. In this case, pk is extended to pk,n, 
denoting the power splitting ratio on subcarrier n at user k. 
Thus, we consider the following optimization problem as 

K N 

cEE Pk,nPnhk,n (20) 

k—1n—1 
K 

'^Xk,n<^yn (21) 

k^l 

a:fc.n € {0,1} ,Vfc,n (22) 

0 </Ofc,„ < l,Vfc, n (23) 

N 

E > C'fc, Vfc, (24) 

n—1 

where rk,n given in dTJ is replaced by 

1 f-I I ~ Pk,n)Pnhk,n\ 

rk,n = log2 ( IH- ^ - I . (25) 

The optimal solution to this problem is given in Appendix 
A. As we have discussed in Section II, (P-UB) is hard 
implemented in currently practical receiver circuits. However, 
it gives a performance upper bound for the comparison purpose 
in simulation. Moreover, since (P-UB) can be directly decom¬ 
posed into several subproblems (details in Appendix A), its 
solutions also provide useful insights to design the step-wise 
algorithm due to its low complexity. 

This step-wise algorithm is executed by two stages. The 
first stage is to seek the optimal subcarrier allocation policy 
X* and the second stage is to find the optimal power splitting 
ratio p*. The two stages are separable instead of correlative 
as in the proposed iterative algorithm. 

The main idea of this algorithm is to first obtain the 
optimal subcarrier allocation variables X* by solving (P-UB), 
then select power splitting ratio p by the bisection search 
individually. We first deduce the following theorem. 

Theorem 1: The optimal subcarrier allocation variables X* 
for (P-UB) is also feasible for (P-PA), given the same inputs 
(same channel conditions and secrecy rate requirements). 

Proof: Please see Appendix B. ■ 

As a result, the proposed step-wise algorithm is feasible as 
long as the optimal algorithm for (P-UB) is feasible. 

Moreover, for each user k, rf. ^ is monotonically decreasing 
in Pk- Therefore, p* can be obtained by the bisection search. 


(P — UB) : max 

{^.p} 

s.t. 
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Algorithm 2 Proposed Step-wise Algorithm for (P-PA) 

1: Obtain X* by solving (P-UB) according to Algorithm 3, 
given the same inputs (the same channel conditions and 
secrecy rate requirements). 

2: for Each user k do 

3: initialize = 1 and = 0. 

4: repeat 

5: Setpfc = i(p^^+p^^). 

6: Compute „ according to O- 

7: if r| > Ck then 

8: Set p^-® = Pfc. 

9: else 

10: Set p^® = Pfc. 

11: end if 

12: until |r^, - Cfcl < eCk- 

13: end for 



The above algorithm is summarized in Algorithm 2. Ac¬ 
cording to Appendix A, the complexity of solving (P-UB) is 
given by 0{K‘^N). In addition, pfc is obtained individually 
by the bisection search, whose complexity is 0{K). Conse¬ 
quently, the complexity of the proposed step-wise algorithm 
is N) which is much lower than that of the proposed 

iterative algorithm. 


IV. Case of Eavesdropper’s Partial CSI 

In this section, we consider a more practical case where only 
statistical CSI of the eavesdropper is known at the BS, i.e., the 
BS only knows the CSI distribution of the eavesdropper. The 
CSI distribution of the eavesdropper can be acquired as follows 
in practice: assuming that the eavesdropper and legitimate 
users are randomly located in the cell (i.e., the eavesdropper 
and legitimate users follow the same distribution), then the BS 
knows the eavesdroppers CSI distribution if the BS obtains the 
legitimate users CSI distribution. Note that this assumption is 
widely used in the literature and more practical than 

the previous case of full CSI of the eavesdropper. 

The performance metric is the ergodic secrecy rate given by 




(26) 


= i log2 1 + 


(1 - Pk)Pnhk,r. 


— Efc 


1 

In 2 


logs 


L 

V 0 - 

l/7e,r 

‘El 


\le,n 


where = ^]E{/3„}, and Ei(a;) = /“ ^dt. 

Note that in our paper pn is assumed to be hxed, the 
eavesdropper information rate re „ in both full and statistical 
CSI cases is independent of pk- Therefore, the analysis of 
both cases is similar. That is, with ([U substituted by (|26]) . the 
problem (P-PA) for the case of eavesdropper’s statistical CSI 
can be solved by Algorithm 1 or 2. 


V. Numerical Results 

In this section, we present the numerical results to evaluate 
the performance of the proposed algorithms. In the simulation 


Fig. 1. Achievable Esum. versus C at total transmit power of 30 dBm. 


setup, we consider an OEDMA network with TV = 128 and 
K = 8 mobile users who are located in a cell of 10 m with 
distance to the BS randomly distributed. The eavesdropper is 
placed exact 10 m away from the BS station. The small- 
scale fading is modeled as i.i.d Rayleigh fading over all 
subcarriers. In addition, the power is uniformly allocated on 
each subcarrier, i.e., p„ = Pt/N, where Pt is the total transmit 
power of the BS. Let TTgum denote the sum power harvested 
by all users. Eor all energy receivers in users’ terminals, it 
is assumed that ( = 0.4. The minimum secrecy rate C is 
assumed to be the same for the all users, i.e., Ck = C,yk. 
Eor the information receivers in users’ terminals, the noise 
power is assumed to be tr^ = —30 dBm. 

Eor performance comparison, we also introduce two 
schemes in simulation as benchmarks. Eor the hrst scheme, 
denoted as hxed power splitting (EPS), power splitting ratio 
Pk = 0.5, Vk, is hxed for complexity reduction and X* is ob¬ 
tained according to (fT4l i. Eor the second scheme, the subcarrier 
assignment is hxed (ESA), while each pk is optimized by the 
bisection search according to Algorithm 2. Specihcally, each 
4 subcarrier is randomly allocated to one user and then we use 
the bisection method to hnd p* achieving all users’ required 
secrecy rate. 

We hrst illustrate the achievable harvested power at different 
required secrecy rate C with total transmit power Pt = 30 
dBm for full CSI case in Eig. [T] It is hrst observed that 
for all schemes, Egum decreases with the increase of secrecy 
rate requirement C. In addition, iTsum falls sharply to zero 
at (7 = 3.51 bit/ OEDM symbol. As we have discussed in 
Section III-B, the optimal X* for (P-UB) can achieve the 
same secrecy rate for (P-PA). Therefore, for both the step¬ 
wise algorithm and the upper bound, Egum falls to zero at the 
same C, where the maximal secrecy rate of both schemes is 
achieved. It is observed that according to the performance of 
the upper bound and the iterative algorithm, applying the same 
power splitting ratio at each user only incurs a little loss in 
terms of the sum harvested power. Moreover, the proposed 
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Fig. 2. Power consumption of information receivers versus C at total transmit 
power of 30 dBm. 

step-wise algorithm incurs less than 35% average loss in 
^'sum compared to the iterative algorithm. Now comparing 
two proposed algorithms with the benckmarking schemes, 
both of them show great advantage over FPS and FSA. In 
addition, the maximal achievable secrecy rate C of the FPS 
and FSA is achieved at around C = 0.45 and 1.5 bit/OFDM 
symbol, respectively, which is much smaller than that of the 
two proposed algorithms. 

Fig. |2] demonstrates power consumption of information 
receivers (the sum received power used to satisfy the required 
secrecy rates) versus different C for full CSI case. We can 
observe that with the increase of the required secrecy rate 
C, more power should be split into the information receivers 
for all schemes. Moreover, the proposed step-wise algorithm 
merely consumes a little more power than the upper bound 
and the iterative algorithm. In addition, the iterative algorithm 
performs close to the upper bound. At last, both proposed 
algorithms consume much less power than FSA and FPS. 

We then illustrate power consumption of information re¬ 
ceivers (the sum received power used to satisfy the required 
secrecy rates) versus the total transmit power Pt for full CSI 
case in Fig. [3] It is first observed that with the increase of 
the transmit power Pt, more power should be split into the 
information receivers for all schemes. This is because with the 
increase of the transmit power, the achievable information rate 
of the eavesdropper ^ also increases, therefore more power 
should also be split into the information receivers to guarantee 
the same secrecy rate. Moreover, the two proposed algorithms 
consume much less power than FPS and FSA schemes. 

In Fig. m we illustrate the relation between iSsum and 
total transmit power Pt with different knowledge of the 
eavesdropper’s CSI (full CSI and statistical CSI) at C = 0.5 
bit/OFDM symbol. First, it is observed that the upper bound 
has the best performance in terms of the sum harvested power. 
In addition, all the schemes are only achievable when Pt > 20 
dBm and perform very close to each other. Moreover, the 



20 22 24 26 28 30 32 34 36 

Transmit Power (dBm) 


Fig. 3. Power consumption of information receivers versus Pt at C = 0.5 
bit/OFDM symbol. 



Fig. 4. Achievable Psum versus Pt at C = 0.5 bit/OFDM symbol. 

upper bound and two proposed algorithms perform better 
with the full eavesdropper’s CSI than only with the statistical 
CSI. However, with the increase of the transmit power, each 
scheme with the statistical CSI performs close to that with full 
CSI, which coincides with the result in 13^ that additional 
channel information of the eavesdropper can hardly provide 
any secrecy rate gain in high SNR region. It is also observed 
that the iterative algorithm performs close to the upper bound. 

We finally demonstrate the relation between 6 and Pt in 
Fig. |5l where <5 is denoted as the ratio of E’sum of one specific 
scheme to that of the upper bound. As it is observed, S for 
proposed iterative algorithm and step-wise algorithm increases 
with the increase of Pt, indicating that the two proposed 
algorithms perform closer to the upper bound with increase 
of the transmit power Pt. When Pt = 37.5 dBm, S of the 
step-wise algorithm achieves 67%, and the iterative algorithm 
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Fig. 5. (5 versus Pt al C = 0.5 bit/OFDM symbol. 


Algorithm 3 Dual-Based Method Algorithm for (P-UB) 

1: initialize A. 

2: repeat 

3: Compute pk,n according to (l36l l and (l37l i. and then ^ 

according to O or (l26l l by replacing pu with pk,n, for 
all k and n. 

4: Solve Xk,n according to dTSl) for all k and n. 

5: Update A via ( |39] | according to the ellipsoid method. 

6: until A converge. 


where A = [Ai, A 2 ,..., is the vector of dual variables. 

The Lagrangian dual function can be obtained as 

where 'R-{X) donate all sets of p for given X that satisfy 
0 < Pk,n < 1 when Xk,n = 1 and pk^n = 1 when Xk,n = 0. 
We can thus obtain the dual problem as 


can reach a 6 = 79.5%. On the other hand, FPS shows no 
improvement in terms of <5 as Pt increases and it always 
performs the worst among all schemes. 

To conclude the discussion on the above results, the pro¬ 
posed iterative and step-wise schemes greatly outperform FSA 
and FPS. Specifically, both carefully coordinating subcarrier 
allocation and selecting power splitting ratio with the adap¬ 
tation to the channel conditions are insignificant improving 
the system performance. Furthermore, while the iterative al¬ 
gorithm performs very close to the upper bound, the step¬ 
wise algorithm also provides favorable performance, greatly 
reducing the complexity. 

VI. Conclusions 

This study investigated the joint subcarrier allocation pol¬ 
icy and power splitting ratio selection for downlink secure 
OFDMA-based SWIPT broadband networks. We formulated 
the problem to maximize the sum harvested power while satis¬ 
fying the secrecy rate requirements of all users. We studied the 
performance upper bound and propose two efficient algorithms 
to tackle the non-convex problems. Numerical results showed 
that the proposed iterative algorithm performs close to the 
upper bound and the proposed step-wise algorithm provides 
a good tradeoff between complexity and performance. 


min q (A) . 


{29) 


The dual function g (A) can be decomposed into N subprob¬ 
lems which can be solved independently. Each subproblem is 
obtained as 


K 


max Lji{Xji^ Pn^ CPn ^ ^ Pk^n^k,n nJ 

X„&T.,p„Gn{X) 

k—1 

(30) 

and we can rewrite the L in (l27l i as 


N 


K 


7 / — ^ ^ Ln ^ ^ 2^kCk- 

n=l fc=l 


(31) 


A. Optimal Power Splitting Ratio 

We first seek for the optimal power splitting ratio of 
each subcarrier. According the Karush-Kuhn-Tucker (KKT) 
conditions ll32l . we have 
1) When rk,n > J'e.n 


dLn 

{^Pk,n 


K 

— ^ ^ CPn^k.n 
k=l 


2^khk^nPn 

ln2 (1 - pk,n) + 


= 0 . 


(32) 


Appendix A 

Optimal Solution for Performance Upper Bound 

We can derive the Lagrangian function for (P-UB) as 
follows: 

L(X,p,A) 

K N K / N 

-EE CPk,nPn^k,n E Ea. E ^k,n'l'k,n ^k 

k—1n—1 k—1 \n—1 

K N K 

=EE (^CPk,nPnhk,n ~t“ ^k,n^k^]^^ri} ^ ^ ^kCf^^ (2V) 

k—1 n—1 k—1 


The optimal solution pk^n can be readily given by 

[ Cln2p„x;f=i^fc.n 

(33) 

where = max{min{*, 6}, a}. 

2) When r^^ri < re,n 


dLn 

dpk,n 


K 

^ ^ CPn^k.n ^ 0- 
k^l 


( 34 ) 



























The optimal in this case can be obtained as 


Pk,n — 1- 


(35) 


Combining the above two scenarios, the optimal solution p’^ ^ 
is summarized as 


Pk,n 


Pk,n^ if 'kk,niPk,n^ ^ e 
1, otherwise, 


(36) 


where 


Pk,n — 


1 - 




n 1 


kf^k,n 


C In 2pn J2k=l ‘^hk^nPn 


■ (37) 


J 0 


B. Optimal Subcarrier Assignment 

Next, substituting p^. ^ into Ln{Xn, Pn), the optimal sub¬ 
carrier assignment policy is given by (the details are easy and 
omitted here). 

^ j 1, if fc = fc* = argmaxfe'Hfc,™ 
x,.„ = < (38) 

’ |0, otherwise, 

where = CPn Ef=l Pk.nhk,n + XkT^n- 

C. Subgradient updating 

As stated in the dual problem is always convex and 
can be solved by using subgradient method. Dual variable A 
can be updated as follow 


On the other hand, for each user k, the secrecy rate 
for (P-PA) is given by 

N 

rlub = 

n—1 

(41) 

which reaches its maximum when pk = 0, and reaches its 
minimum when = 1. 

Thus, for the given set of X, we obtain that 


logs 1 + 


(1 Pk') Pn^k^r^ 


- Te 


i,{pk,n = 0) 

Pk,n 


N 

— ^ ^ Xk,n 


1 ( 1 I Pnhk,r. 

I0g2 ( 1 + —^ 


= maxrf = rf (pfc = 0). 


Pk 


(42) 


In another word, the maximal secrecy rates of both case 
equal, given set of X. 

Furthermore, for feasible solution X* and {p^^} for (P- 
UB), we have 


0 < {plj) < Tl^,{X\ {p,,„ = 0}) 

= r^,pJX*,pfc=0). (43) 

Since is a continuous function with respect to p^ and 
monotonically decreasing in p^, there always exists a certain 
Pk € [0,1] that satisfies 

0 < rl,AX*,pk) = Tl^,{X\ {plj) < rl^,{X\pk = 0). 

(44) 




A1*^ + Q-k 


t N 

Ck ^ ^ Xk^', 


,yk. (39) 


The complexity of this dual based algorithm is analyzed as 
follows. For each subcarrier, 0{K) computations are needed. 
Since the calculation is independent at each subcarrier, the 
complexity if 0{KN) for each iteration. Last, the complexity 
of subgradient updates is polynomial in K ll^ . Hence, the 
overall complexity of subgradient method is N). 

Finally, we present the whole algorithm in Algorithm 3. 


Appendix B 
Proof for Theorem 1 


In this appendix, we will prove that the optimal subcarriers 
assignment X* for (P-UB) is also feasible for (P-PA), given 
the same inputs (the same channel conditions and secrecy rate 
requirements). 

On one hand, the secrecy rate for user k for (P-UB) is given 
by 


N 


' k,ub 


= E 


^k,n 


n—1 


l0g2 1 + 


(1 Pk,n) Pn^k.T 


1 + 


(40) 

where reis the information rate of the eavesdropper on 
subcarrier n given in ([T]i (full CSI) or (l26l l (statistical CSI). 

For each user k, with the fixed feasible X reaches its 
maximum when pk,n = 0 for all n that satisfy Xk,n = 1- 


The proof is thus completed. 
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